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Abstract 

Motivated by Hcilmann and Licb's work [H], we discuss energy level crossings for the 
one-dimensional Hubbard model through the Bethe ansatz, constructing explicitly the 
degenerate eigenstates at the crossing points. After showing the existence of solutions 
for the Lieb-Wu equations with one-down spin, we solve them numerically and construct 
Bethe ansatz eigenstates. We thus verify all the level crossings in the spectral flows 
observed by the numerical diagonalization method with one down-spin. For each of the 
solutions we obtain its energy spectral flow along the interaction parameter U. Then, 
we observe that some of the energy level crossings can not be explained in terms of U- 
independent symmetries. Dynamical symmetries of the Hubbard model are fundamental 
for identifying each of the spectral lines at the level crossing points. We show that the 
Bethe ansatz eigenstates which degenerate at the points have distinct sets of eigenvalues 
of the higher conserved operators. We also show a twofold permanent degeneracy in terms 
of the Bethe ansatz wavefunction. 

1 Introduction 

Degeneracies in the energy spectra of quantum systems have close relationships with their 

symmetries. Actually, von Neumann and Wigner showed that degeneracies are more likely to 

occur for the systems with one or more symmetries than those without symmetries |2U1 113j . 

To be more precise, if one assumes that a quantum system is given by a real Hamiltonian 

matrix whose elements are expressed by independent parameters, in the case of no symmetry, 

two parameters happen to take some prescribed values in order to bring two of the eigenvalues 
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into coincidence. Their theory reminds us of the "non-crossing rule" in quantum chemistry, 
which states that energy levels of orbitals of the same symmetry can never cross each other 
along a reaction parameter. However, von Neumann-Wigner's theorem does not give a proof 
for the non-crossing rule. It is possible that degeneracies appear in the systems without sym- 
metries. Such degeneracies are referred to as accidental degeneracies. In fact some examples 
of accidental degeneracies are numerically observed in molecules j21| or triangular quantum 
billiards p. 

The one-dimensional Hubbard model is one of the most significant models in condensed 
matter physics. The model also attracts a great interest of mathematical physicists due to 
its Bethe-ansatz solvability jlUl [Hj . Heilmann and Lieb numerically investigated energy spec- 
tral flows along the interaction parameter U for the system on a periodic 6-site chain and 
found many level crossings which can not be accounted for by the known symmetries such 
as translation, SO(4) and particle-hole symmetries [Sl|2j. They concluded that, if one takes 
into account only [/-independent symmetries, the level crossings are accidental degeneracies, 
that is, a counter example of the non-crossing rule. Recently, Yuzbashyan, Altshuler and 
Shastry have suggested that the origin of Heilmann-Lieb's level crossings should be dynam- 
ical symmetries in the Hubbard model |2H|. Here the dynamical symmetries are given by 
parameter-dependent operators, which are often called higher conserved operators in associa- 
tion with conserved quantities in classical integrable systems. The dynamical symmetries for 
the Hubbard model are constructed in ^| El HI l22l • Yuzbashyan, Altshuler and Shastry 
numerically showed that crossings in the spectral flows of the first three conserved operators 
never occur at the same value of U. The dynamical symmetries depend on the parame- 
ter U, and they are not considered as symmetries in the von Neumann-Wigner's theorem. 
Heilmann-Lieb's level crossings are still considered to be accidental degeneracies. 

In the framework of the Bethe ansatz, we discuss in the paper energy level crossings for 
the one-dimensional Hubbard model. The Bethe ansatz method provides information on the 
eigenstates that can not be easily obtained in the direct diagonalization of the Hamiltonian 
matrix. We set up the following problems: 

i) When two energy eigenvalues approach in numerical data as one parameter is varied, one 
may draw two alternative spectral flows, a level crossing or level repulsion |S]. To ensure 
that genuine energy level crossings have happened, we must investigate the change of 
each eigenstate along the parameter U. 

ii) Do the eigenstates have distinct dynamical symmetries at the level crossing points? In 
order to solve the problem we assign each of the eigenstates the eigenvalues of the higher 
conserved operators. 
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It is indeed not easy to investigate these problems. For the triangular quantum billiards with 
two parameters, Berry and Wilkinson investigated the behaviour of eigenstates along a circuit 
of the crossing point in the parameter space so that they could verify the existence of genuine 
energy level crossings 0. For general quantum systems, it is hard in practice to construct 
quantum many-body eigenstates. The task is also not easy even for the systems that can 
be treated by the Bethe ansatz. In fact, it is nontrivial to obtain numerical solutions to the 
Bethe ansatz equations for a finite size system. For the sector of one down-spin, however, it 
is practically possible to solve the Bethe ansatz equations numerically. 

In the present paper, we prove the existence of solutions to the Lieb-Wu equations and 
then numerically solve them. Here we generalize the method of Ref. By using the 
numerical solutions, we analyse the behaviour of each of the degenerate eigenstates along the 
parameter U . We thus find several genuine level crossings of energy spectral flows with the 
same [/-independent symmetries. The merit of this method is that the numerical solutions 
for the Lieb-Wu equations provide not only the eigenvalues of the higher conserved operators 
but also the one-to-one correspondence between their spectral lines and eigenstates. As a 
consequence, we observe that all the common eigenspaces of the first three higher conserved 
operators are one-dimensional in the subspaces with the same [/-independent symmetries. 
Furthermore, when there are [/-independent degeneracies (permanent degeneracies) in the 
spectral flow, the one-to-one correspondence plays an essential role in assigning to each of 
the degenerate eigenstates its eigenvalues of the dynamical symmetries at the energy level 
crossing points. By using the explicit form of the Bethe ansatz wavefunctions 23^^, we can 
derive the permanent degeneracies. We remark that the existence of level crossings also gives 
a necessary condition for algebraic independence of the three higher conserved operators in 
the subspaces. 

This paper is organized as follows: first, in Section]^ we summarize the known results for 
the Hubbard model and prepare for the following sections. We list the known symmetries of 
the Hubbard model in 12.11 and review the Bethe ansatz method in 12.21 In 12.31 we prove the 
existence of solutions for the Lieb-Wu equations with one down-spin following the approach 
developed in Next, in Section |31 we analyse degeneracies in the energy spectrum of the 
Hubbard model by using both direct diagonalization of the Hamiltonian matrix and the Bethe 
ansatz method. In 13.11 we describe twofold permanent degeneracies arising from a reflection 
symmetry of the lattice. In 13.21 and 13.31 we investigate energy level crossings for the systems 
with two or three up-spins and one down-spin on a periodic 6-site chain. We find some 
genuine energy level crossings and see that the degenerate eigenstates can be classified by 
the eigenvalues of the higher conserved operators. The final section is devoted to concluding 
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remarks. 



2 Symmetries and Bethe ansatz method 

We introduce the Hubbard model on a one-dimensional periodic L-site chain. Let cj^ and 
Cis,{i € 'LjUL^s S {T) i}) be the creation and annihilation operators of electrons satisfying 
{ci^5,Cjt} = {cj^,c^^} = and {cj^jC^^} = Sij6st, and define the number operators by riis := 
clgCis- We consider the Fock space of electrons with the vacuum state |0). The one-dimensional 
Hubbard model is described by the following Hamiltonian acting on the Fock space: 

^ = - E E (4ci+l,. + cUl,sCis) + UY^ ("it - 2) ("^i - 2) ' (^-^^ 
i=l s=],l i=l 

where C/(€ M>o) is the interaction parameter. We consider the system for L even throughout 
the article. 

2.1 [/-independent and dynamical symmetries 

In general, the symmetries of a quantum system are expressed by operators which commute 
with its Hamiltonian. They are classified into two families in the case of the Hubbard model; 
one is independent of U and another depends on U. We list some of these symmetries. First 
we consider the [/-independent symmetries 0. Define 

- - n n - ■■= n ^if^s ■ ■ ■ ^ii,.' (2-2) 

where pj^^^ := 1 — (c|^ — CjJ{cis — cjs), {i 7^ j) are permutation operators. The operators a 
and T correspond to reflection and translation symmetries of the lattice, respectively. They 
satisfy = 1. It is clear that they commute with the Hamiltonian (|2.1|) . Another 

[/-independent symmetry is the SO {A) symmetry [251 126j . Define 

^ L L 

Sz ■■= ^^{rii^ -nil), S+ ■.= ^cl^Cii, := {S+)\ 

1=1 1=1 
^ L L 

Vz ■■= -^{l-nn -nil), V+ ■=J2^~y^'^^^^' ri--={ri+)^- (2-3) 

1=1 i=l 

Both sets of operators {Sz,S±} and {r]z,r]±} give representations of the algebra su(2) in the 
Fock space. They all commute with the Hamiltonian 1)2. which leads to the symmetry of 
type su(2) © su(2) ~ so(4). Furthermore it is known that this so(4) symmetry lifts to the 
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50(4) group symmetry. For the later discussion, we define Casimir operators, 

52 := i + S^S+) + Si r/2 := 1 (r?+r?_ + + tjI 

which commute with all the operators in (|2.3j) . 

Next we introduce the dynamical symmetries given by the [/-dependent operators. These 
[/-dependent operators themselves are also commutative and are called conserved operators 
in association with conserved quantities in classical integrable systems. In |15| 17]. the first 
three conserved operators are explicitly given by 

h=H, 

L 

^2 = - 1^ X] (4ci+2,.-cl+2,sCis) 
L 

+ J2 (4ci+l,s-4+l,sCis)(^i-^ + "'i+l -s-l)' 

L 

^3 = - XI (4ci+3,. + 4+3,sCis) 

- (n..-i) (nm,-.-^) -li^^s-l) {n.,~s-l)) 

In order to obtain higher conserved operators systematically, the transfer matrix approach 
similar to that of the XXZ Heisenberg spin chain is developed |141 1151 VFR lllj . The 50(4) 
symmetry of such higher conserved operators is also verified in the framework of the transfer 
matrix approach |16j . 

The operators {r, S^, 5^, 77^, 772, /„} give a commutative set of operators including the 
Hamiltonian. Hence they can be diagonalized simultaneously. 

2.2 Bet he ansatz method 

The Bethe ansatz method was applied to the Hubbard model in 10 . Here we give only 
the result. Let be the number of electrons and M that of down-spins. We may assume 
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2M < N < L due to particle-hole and spin reversal symmetries in the system. Let k = 
{ki\i = 1, 2, . . . , A^}, (Re(fct) € M/27rM) denote a set of wavenumbers of electrons and A = 
{Aala = 1,2,...,M} that of rapidities of down-spins. Given a set of spin configuration 
s = {si\i = 1, 2, . . . , N} with N — M up-spins and M down-spins, the Bethe state with k and 
A has the following form: 

\k,\:s) = ^il^k,x{x]s)cl^^,^cl^^,^ - ■ ■ cl^^,J{)). (2.4) 
The coefficients 'ipk,\ix;s) in ()2.4|) are expressed [15] as 

N 

i'kAx; s) = ^ sign(PQ)(/7fe^,A(sQ)exp (^i^fcp(i)XQ(j)^, 

Pe6]v i=l 

Ej—r- \fi(ce) — ^R(f3) — lf^/2 y-r 
11 — T ^ \\Fkp{XR{^),y^), 

p , 1 Xa-smkp(j^-iU/A 

A„-sinA;p(j^)+iC//4 Ac,-sinyfcp(j)+i?7/4' 

where we have denoted by Q one of the shortest elements in the symmetric group Sat on 
{1,2,..., A^} such that 1 < < xq{2) < • • • < < L, and by the position of the 

7th down-spin in sq = {sq(i), sq(2)) • • • > ^(^{Ar)}. The Bethe states (|2.4() give eigenstates of 
the Hamiltonian 1)2. ip if {/cj, Aq,} satisfy the following equations: 

g^^fc^L ^ TT A/3-sinfc^-i[//4 
A/3-sinA:, + i[//4' 

Aq, - sinkj - \U/A _ -j-j- Aq, - A;3 - \U/2 
.^^ Aq - sin ki + i?7/4 " ^ Aq - A;3 + iU/2 ' ^ ' ^ 

which are coupled nonlinear equations called Lieb-Wu equations. The Lieb~Wu equations 
have not been solved analytically. But it predicts some important results on thermodynamic 
properties of the system through Takahashi's string hypothesis |17| I18 | 1^1 119j. 

The Bethe states (|2.4j) are not only eigenstates of the Hamiltonian ()2.1j) but also those 
of the translation operator r and the higher conserved operators I2 and I3. By using the 
solutions {A;i,AQ,} of the Lieb~Wu equations (|2.5|) . the eigenvalues of r and {/„} are written 
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as 

L ^ 

T\k,\;s) = e'T-P|/c, A; s), P = ^(X/^*) (™odL), 

i=l 

In\k,\]s) = En\k,\]s), (n = l,2,3), 

^ 1 
Ei = E = -2^cos ki + -U{L-2N), 

i=l 

N 

E2 = -2^{ sm{2ki) + U sin ki) , 

E3 = -^(2cos{3ki) + 3u(cos{2ki) - +U'^cosk^^ - -UL. (2.6) 

i=l 

The P appearing in the eigenvalues of r indicates the total momentum of the system. (There 
should be no confusion with the use of P in the coefficients of the Bethe states where P 
denotes an element in &n-) 
We immediately find 

S,\k,X;s) = ^{N-2M)\k,X;s), r],\k, X; s) = ^{L - N)\k, \; s). (2.7) 

It is shown in [Hj that each Bethe state (|2.4|) with a regular solution for the Lieb-Wu equa- 
tions (|2.5j) corresponds to the highest weight vector of a highest weight representation of 
so(4), i.e., S-^-\k, A; s) = ri^\k, A; s) = 0. Then we find 

S^\k,X;s) = S{S + l)\k,X;s), ri'^\k, X; s) = r]{7] + l)\k, X; s), (2.8) 

with S = {N — 2M)/2 and r] = (L — N)/2. Hence, by applying the lowering operators 
(5_)", (0 < n < iV - 2M) and (r?„)™, (0 < m < L - iV) to the Bethe states ^T^, we also 
have eigenstates of the Hubbard Hamiltonian 1)2. 

\k,X;s;n,m) := (5_)"(r7_)"*|A;, A; s), (2.9) 

which have the same eigenvalues {En} for the operators {In} as those of \k, A; s). The dimen- 
sion of the representation with the highest weight vector \k, A; s) is {N — 2M -|- 1)(L — -|- 1). 
By using Takahashi's string hypothesis for the Bethe states 1)2. 4|) together with the 50 (4) 
symmetry, their combinatorial completeness is proved in 0. 

2.3 Lieb— Wu equations with one down-spin 

We try to find regular solutions of the Lieb-Wu equations ()2.5() in the case when the system 
has only one down-spin following the discussion in [21. In this case, the string hypothesis ^7] 
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predicts that two types of solutions exist: one is the solution with only real wavenumbers 
{ki} and another includes two complex wavenumbers. First we consider the real wavenumber 
solutions. For M = 1, the Lieb-Wu equations 1)2. 5|) reduce to 



_ A-sin/ci -1C//4 



X-sinki + iU/A' ^ ' ' . ^ 



These are equivalent to the following equations: 



U k-L ^ 
sin fej — A = — cot ( — ^ j , exp (1 kjLj = 1. (2-11) 



i=l 

We consider the real solutions for the first equation 

sing-A = ^cot (^). (2.12) 
In the interval < q < 27r, its right hand side has L branches 

|(^4)<'<l(^4)' (2.3) 

If we seek a solution q of (|2.12|) in one of the branches (|2.13jl , the solution is unique under the 
condition UL > 8 and can be written as a function of A, i.e., q = qi{X). Given a distinct 
set = 1,2,..., N} C = 1, 2, . . . , L} of the branches, the second equation in 1)2. 12() 

is satisfied when ^ J2iLi QiiW ^ ^- The behaviour of the solution q = qi{X) tells us that 

L ^ ^ 1 

i=l i=l 



Hence there exist — 1 values of A which give the following integer values for T^i=i Qiii^)'- 

N 

2. 



m e 



N 



1=1 



i = l,2,--- ,iV-l 



Note that such {A} and integers {m\ are in one-to-one correspondence due to -^^^ > 0. It 
is straightforward that {ki = q£-{X),X} characterized by the indices {£i,m} give (^)(-^ — 1) 
solutions of the equations 1)2. 11() . 

Next we consider the fc-A-string solutions. We assume the forms of solutions as 

ki G M/27rM, (i = 1, 2, . . . , - 2), kN-i = C-i?, kN = C+iC, 
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where < ^ < 27r and > 0. Note that kjy-i and k]\f form a complex conjugate pair which is 
referred to as fc-A-two-string. Then the first equations in (|2.inj) are rewritten as the foHowing 
equations with real variables: 

U /kL\ 

sin/ci - A = — cot (^), = l,2,...,iV-2), (2.15a) 



4 V 2 

smCcosh^ - A = J , "f^"^^ (2-15^) 
4 cosh(4L) — cos(C-/^) 

, ^ U sinh(fL) 

'=°"^''°''^ = -4 coshKL)-cosKL) - 

On the other hand, the second equation in (|2.1()|) is equivalent to the following condition: 

\^ki + 2Q = —m, (mG{0, 1,...,A^L-1}). (2.16) 

i=l 

In the same way as the previous case, if we consider a solution of each equation H2.15a() 
in one of the branches (|2.1l'{j) . it can be written as a function of A. Given a distinct set 
= 1,2, . . . , N — 2} of indices specifying the branches (|2.13|) . we express the solutions 
of (I2.15al) as ki = qiX>), {i = 1,2, . . . , N - 2). Then, from the relation (ITTKl) . the C is also 
written as a function of A, 

^ N-2 

( = ((X):=^m--Y,Qi.W, (2.17) 

i=l 

for fixed {£{} and m. For an illustration, we consider (|2.15b|) and (|2.15c|) in the case N = 2. 
Since C, does not depend on A in the case, the equations (|2.15bl) and (|2.15c|) decouple into 

A = sin(^m)coshe, sinhg = - ^ J ^'HO, (2-18) 
\L J 4cos(.^m,j 



where 



(2),^^ sinh(^L) J tanh(^L/2) for m odd. 



cosh(^L) - (-I)*" [coth(CL/2) for m even. 

One finds from graphical discussion 3 that, if the condition C/L > 8 is satisfied, the second 
equation determines an unique ^(> 0) for ^ < < We denote it as ^ , {m G = 
1, 2, • • • , L— 1}). The first equation in 1)2. 18() immediately gives A with the . Let us consider 
the case N > 2. By inserting (|2.17|1 into (|2.15c|l . we have 

sinhe = -— ^^/(e), /(O := — — — — . (2.19) 

4cos(C(A)) cosh(eL) - (-)*" cos (f ^^(A)) 

Since /(^) > for > in the similar to the case N = 2, this determines an unique ^ as a 
function of A if and only if f < ({X) = f m -^Y^i 9^, (A) < By using 1)2. 14|) . it is sufficient 



9 



to have an unique solution for ()2.19() that the integer m satisfies 

N-2 , ^ N-2 



i=l i=l 



that is, 



I 1=1 



1,2,...,L-N + 1}. (2.20) 



(2) 

Note that hmA-+±oo = C \- a, which is weh-defined for the above m. We expect 
that, for the L — N + 1 values of m allowed in H2.20() . the equation H2.15b|) with ^(A) and C{X) 



A = sin(fm-iX:.'?^.(A))cosh (e(A)) 



4cos(f E,,%(A))-(-rcosh(e(A)L) 
= :g{qUX),aX)), (2.21) 

determines A. In fact, since qi^iX) and ^(A) are continuous functions with respect to A and 

^lim^,(,,(A),e(A))=,(f(^.±i),ej:^,.,.)), 

5 is a continuous and finite function with respect to A. Hence there exists a solution A in 
(|2.21|) . As a consequence we have (^^g)!-^ — + 1) solutions corresponding to the indices 
{ii,m}. 

Let us see the relation between the string hypothesis ^7j and our results. Let {£{} reduce 
modulo L into the interval [—■§,-§) and express them as {/^*^}. Putting J = m — J2iLi 
the real wavenumber solutions and J' = L — m + fc-A-two-string solution, we 

have 

jWez + i, -:^<i»<:^, 
2' 2 - 2' 

A^ 1 A^ 1 

JGZ + -, \J\<-{N-2), J'GZ+y, \J'\<-{L-N). (2.22) 

One sees that the indices {I^*^; J; J'} characterizing the regular solutions of the Lieb-Wu 
equations with M = 1 are nothing but those appearing in the string hypothesis J7]. Thus we 
have shown for the Lieb-Wu equations 1)2. 5|) with M = 1 that there exist the same number 
of solutions as those predicted by the string hypothesis. In the next section, we numerically 
calculate the solutions {ki, A|z = 1, 2, . . . , A^} and {/cj, C, ^, A|i = 1, 2, . . . , A^ — 2} for L = 6. 

3 Energy level crossings 

Let us review von Neumann-Wigner's discussion on the spectra of quantum systems. We 
assume that a Hamiltonian is described by a finite-dimensional real symmetric matrix whose 
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elements are regarded as random independent parameters. If the system has no symmetry, 
we call its spectra a "pure sequence". If there exist some symmetries in the system, its 
spectra is given by a superposition of pure sequences, which we call a "mixed sequence" . Von 
Neumann-Wigner's theorem reads as follows: one must adjust two parameters to bring two of 
the eigenvalues belonging to a pure sequence into coincidence while in a mixed sequence one 
obtain a degeneracy by varying only one parameter |2U1 113j . Hence the degeneracies in pure 
sequences are very unlikely to be found if we choose the values of parameters in an arbitrary 
manner. Such degeneracies in pure sequences are referred to as accidental degeneracies. 

Applying the above discussion, we study the Hubbard model. In varying the parame- 
ter U, the Hamiltonian (|2.1|) gives a flow in the above space of parameters. As we have 
mentioned in the previous section, the system has several symmetries. Since the opera- 
tors {t, S , Sz, rf,r]z} are mutually commutative, they can be simultaneously diagonalized 
by an orthogonal transformation. Through the same orthogonal transformation, the Hamil- 
tonian matrix breaks up into diagonal blocks corresponding to the common eigenspaces of 
{r, S^, 5*2, r/^, Notice that the common eigenspaces are characterized by the set of quan- 
tum numbers {A^, M, P, 5, r/}. Energy eigenvalues from the blocks with different quantum 
numbers may degenerate due to translation and 50(4) symmetries. But all the blocks do not 
give a pure sequence; for example, the blocks with L = N = 2M have particle-hole and spin 
reversal symmetries. However, after considering all the known [/-independent symmetries, 
the spectra also have many degeneracies at special values of [/, i.e., level crossings in the 
spectral flows along the parameter U ;8 . Thus the flow determined by the Hamiltonian (|2.1() 
runs through several special values of parameters that give accidental degeneracies. 

Let us discuss how to numerically determine level crossings, in particular, for the case of 
accidental degeneracies. One notices that the numerical diagonalization of the Hamiltonian is 
not enough to find crossings of energy spectral flows since apparent crossings may be just close 
approaches of two eigenvalues |Hj. To verify the existence of genuine level crossings, we must 
investigate the behaviour of spectral flows for each eigenstate. The Bethe ansatz method is a 
very effective tool to establish this. Here we restrict our investigation to the systems with one 
down-spin where we have shown the existence of solutions for the Lieb-Wu equations 1)2. 5|) in 
Subsection 12.31 We have seen in Subsection 12.11 that the Hubbard model has the dynamical 
symmetries in addition to the [/-independent symmetries. The recent paper |2H1 pointed out 
that the degenerate eigenstates at the accidental degeneracies observed in |H] can be classified 
by the eigenvalues of the higher conserved operators. We verify their assertion at the level of 
eigenstates in the sector of one down-spin. Here we note that the degeneracies discussed in 
[HI 121 128j are in the half-filled Hubbard model with zero magnetization, that is, in the sector 
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of three down-spins. The strategy in the present paper is given by the following: 

i) We give a matrix representation of the Hamiltonian in a certain common eigenspace of 
the operators {r, S'^, 5^, ry^, r/^}, which is referred to as desymmetrization of the corre- 
sponding symmetries. We numerically give its eigenvalues from the direct diagonaliza- 
tion [HI US] for several values of U. 

ii) We give the numerical solutions of the Lieb~Wu equations ()2.5|) in several values of U. 
From the correspondence between energy spectral flows and the Bethe states, we verify 
whether or not genuine energy level crossings exist. 

iii) We also diagonalize other conserved operators {12,13} and see the correspondence be- 
tween their spectral flows and the Bethe states. We analyse the structure of their 
degeneracies. 

We discuss only the systems containing two or three up-spins and one down-spin, which do 
not have particle-hole and spin reversal symmetries. In spite of such restriction, we have some 
non-trivial results on degeneracies. 

3.1 Twofold permanent degeneracies 

The translation and the 50(4) symmetries produce [/-independent degeneracies in energy 
spectral flows. We call [/-independent degeneracies permanent degeneracies. Furthermore, 
after the desymmetrization of these [/-independent symmetries, we often observe another 
twofold permanent degeneracies associated with a reflection symmetry of the lattice [2Sj. In 
fact, we can explain them in terms of the Bethe ansatz wavefunction. Here we note that, due 
to the relation ra = crT~^, the reflection operator a acts only on the eigenspaces of r with 
the eigenvalue 1 or —1, that is, the subspaces of the Fock space with the total momentum 
P = or |. 

Let us investigate the twofold permanent degeneracies due to the reflection operator a 
in the framework of the Bethe ansatz method. Even in the subspaces with P = or l^, the 
Bethe states do not always diagonalize the operator a. Indeed it is easy to verify that, if we 
apply the operator a to one of the eigenvectors \k, X; s;m,n), then its total momentum and 
eigenvalues of {hn} are negated and those of {S'^ , Sz,T]'^,'riz, hn-i} do not change since 

ar = T~^a, ahn + hnf^ = 0, ahn-i = hn-icr, [a, Sz,±] = [cr, i]z,±] = 0. (3.1) 

These facts are verified more directly from the following relation: 

a\k, A; s; n, m) = (— )^'^| — k, —A; s; n, m), (3.2) 
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where —k = {—ki\i = 1, 2, ... , N} and —A = {—Xa\a = 1, 2, ... , M}. Note that, if {ki, X^} 
is a solution for the Lieb-Wu equations (|2.5|1 . then so is {—ki,—Xa}. The relation (|3.2j) 
means that, if the set {k^; Xa} does not coincide with {—ki;—Xa} as a set, we have a twofold 
permanent degeneracy in energy spectral flows. On the other hand, if {ki; Xa} = {—kf, —Xa} 
as a set, then the eigenstates \k, X; s; m, n) have the eigenvalue £"2 = for the second conserved 
operator I2, which follows from the formulas in ()2.6() . It should be noted that, to see the 
existence of such twofold permanent degeneracies, we must solve the Lieb-Wu equations 1)2. 5() . 

3.2 Spectral flows: L = 6,N = 3,M = 1 

We now study the system with two up-spins and one down-spin on benzene {L = 6, N = 
3, M = 1). Note that Sz = and t?^ = | in this case. We consider the subspace characterized 
by the set of quantum numbers {P = 2, S = ^,7] = |}. There is no more [/-independent 
symmetry in the subspace, and we shall actually find one of the simplest nontrivial energy 
level crossings there. We have a matrix representation of the Hamiltonian of 11 dimension. 
We put u := U /{U + A) and numerically diagonalize the Hamiltonian matrix for several values 
of u. FigureHshows the energy spectral flows for < u < 1. Here the vertical line in the figure 
indicates the energy E divided by C/ + 4. We confirm from the numerical data that there is 
no permanent degeneracy. We observe a close approach of two energy levels in 0.5 < u < 0.6, 
which seems to be an energy level crossing. 

Our first purpose is to show that a genuine energy level crossing has been found in Figure 
n In the case of triangular quantum billiards , topological properties of their eigenstates 
were studied to verify the existence of energy level crossings. Here we employ the Bethe 
ansatz method to verify energy level crossings in the level of eigenstates. It is clear that, since 
S = Sz = 12 and = ry^ = |, all the eigenstates in the subspace are the Bethe states charac- 
terized by the indices {/^*^; J; J'} satisfying (|2.22|) . By using the procedure in Subsection l2.3[ 
we numerically give real wavenumber solutions with P = (Y^^i Z^*-* + J) (mod 6) = 2 (mod 6) 
and /c-A-two-string solutions with P = (/^^^ — J') (mod 6) = 2 (mod 6) to calculate the corre- 
sponding energy eigenvalues. The correspondence between energy eigenvalues and the Bethe 
states at u = 0.3 and 0.8 is displayed on Table 1. Here we deal with the Lieb~Wu equa- 
tions only in the case when the condition UL > 8 is satisfied. One sees that the results for 
u —> 1,{U — > 00) read as 

E I for real wavenumber solutions, 

lim = < 

U + 4 I 1 for A;-A-two-string solutions, 

which agrees with those conjectured by the string hypothesis. We remark that the energy 
eigenvalues obtained by the solutions of the Lieb~Wu equations coincide with those obtained 
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by the direct diagonalization of the Hamiltonian matrix within an error of 0(10^^^). Thus, 
by combining the results in Figure ^ with those on Tabled we obtain the behaviour of energy 
spectral flows for each eigenstate. We conclude that, in 0.5 < u < 0.6, there exists an energy 
level crossing between two Bethe states indexed by {— |, |, |; ^; } and ; — |}. 

Next we show that, if we take into account dynamical symmetries, each one-dimensional 
component of the degenerate eigenstates at the energy level crossing point can be distin- 
guished from the other degenerate eigenstates. In the similar way to the energy eigenvalues, 
the eigenvalues of the second conserved operator I2 are obtained by both their direct diago- 
nalization and the solutions of the Lieb-Wu equations, which are displayed in Figure [21 and 
Table 121 (the vertical line in the figure also indicates the eigenvalue E2 divided by 17 + 4). 
Note that there is no permanent degeneracy. We observe that the spectral flows indexed 
by {— and {^;;— |} in Figure [21 never have crossings. Hence the Bethe states 
indexed by {— |, |, |; ^; } and ; — |} belong to different eigenspaces of I2, that is, the two 
Bethe states have a different dynamical symmetry. Thus all the common eigenspaces of the 
operators {H,l2} with the set of quantum numbers {P = 2, S = Sz = ^,ri = rjz = ^} are 
one-dimensional. 

The fact that the common eigenspaces of {H,l2} are one-dimensional becomes clearer by 
investigating the spectral flows of the transfer matrix, i.e., the spectral flow sets of conserved 
operators. Figure [SI shows that, if we consider the eigenvalue sets {£,£2} for two of the 
eigenvectors along the parameter U, there is no crossing while their projection onto the 
U-E plane has a crossing. One notices that the Bethe ansatz method plays an essential 
role in getting such picture. It seems to be rare that the eigenvalues of the transfer matrix 
simultaneously degenerate. We comment that the XXZ Heisenberg spin chain at root of unity 
has such degeneracies in the transfer matrix where one sees symmetries defined only at 
the special values of parameters. 

3.3 Spectral flows: L = 6, = 4, M = 1 

We give an example where the eigenvalues E and E2 degenerate at the same value of U. We 
consider the system with three up-spins and one down-spin on benzene {L = 6, N = 4:, M = 1). 
In this case, we have Sz = ijz = 1. We investigate the degeneracies of energy spectral flows 
in the subspace characterized by the set of quantum numbers {P = 0, S = rj = 1}; hence we 
can expect all the eigenstates therein to be the Bethe states. The subspace is 14 dimensions. 
Figure 01 shows the energy spectral flows obtained by the numerical diagonalization of the 
Hamiltonian matrix for < u < 1. Also listed on Table [SI are the energy eigenvalues given 
by numerical solutions for the Lieb-Wu equations at u = 0.4 and 0.8; we show the numerical 
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Figure 1: Spectral flows E from direct diagonalization for L = 6, A'^ = 3, M = 1, P = 2, S* = ^ 
and ?7 = I ■ 
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Table 1: E obtained by the Bethe ansatz method. 
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Figure 2: Spectral flows E2 from direct diagonalization for L = 6, = 3, M = 1, P = 2, S* = | 
and ?7 = I ■ 
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Table 2: E2 obtained by the Betlie ansatz method. 
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Figure 3: Spectral flow sets (i?, £"2) characterized by {— and {^;;— |} for the 
system with L = 6,N = 3,M = 1, P = 2,S = ^ and ?? = |, and their projections onto the 
U-E and U-E2 planes. 

solutions in Appendix \K\ We observe as u ^ 1, {U 00) that 

E I — i for real wavenumber solutions, 
lim — = < 1 

U—^oo (7+4 I - for A:- A- two-string solutions, 

which are also derived from the string hypothesis .l7j. Here we note that the results from 
the two different methods above coincide within an error of O(10~^^), which give an evidence 
for the validity of the Bethe ansatz method. From Figure |1] and Table IHl we see the corre- 
spondence between the energy spectral flows and the indices {Z*-*^; J; J'} (|2.22j) characterizing 
the Bethe states. The correspondence shows that, if J; J'} 7^ {— I*^*); — J; — J'}, i.e., 
[ki] A} 7^ {—ki; —A} as a set, two Bethe states characterized by these indices are in a perma- 
nent degeneracy. In fact, among the 14 Bethe states, 8 of them are in the twofold permanent 
degeneracies associated with the reflection symmetry. We find that the Bethe states in- 
dexed by {— |, — |, |, |; 0; } and {— ^, ^; ; 0} has an energy level crossing at n = 0.5, {U = 4). 
Although the existence of this energy level crossing is also verified from the characteristic 
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equation of the Hamiltonian matrix, we need the Bethe ansatz method to clarify which Bethe 
states have the crossing. 

We discuss whether or not each of the above degenerate eigenstates have distinct eigen- 
values of the higher conserved operators. In the similar way to the previous case, we consider 
the correspondence between the spectral flows of the second conserved operator I2 and the 
Bethe states, which is displayed in Figure [31 and Tabled As is described in Subsection l3.11 if 
we express the eigenvalue of I2 for the eigenstate \k, A; s) as E2, that for the state a\k, A; s) is 
—E2. In particular, if {I*^*^; J; J'} = {—/'•*•'; — J; — J'}, i.e., {h; A} = {—ki] —A}, the state has 
E2 = 0. FigureEland TableElindeed show that two eigenstates whose energy spectral flows are 
in a permanent degeneracy belong to the different eigenspaces of I2 and there are 6(= 14 — 8) 
states in the eigenspace of I2 with E2 = 0. Thus we may conclude that the dynamical sym- 
metry I2 "accounts for" all the twofold permanent degeneracies in the energy spectral flows 
which are produced by the reflection symmetry. It is a further remarkable fact that both 
the Bethe states indexed by {— |, — |, |, |; 0; } and {— ^,^;;0}, which has an energy level 
crossing, belong to the same eigenspaces of I2 with E2 = 0. Hence, at n = 0.5, {U = 4), the 
two Bethe states can not be distinguished from their eigenvalues of {r, S'^ , SzjT]"^ jij^, H, I2}. 

We have to investigate the spectral flows of higher conserved operators. Displayed in 
Figure IHl and Table [5] are the spectral flows of the third conserved operator obtained 
by its direct diagonalization and its eigenvalues at u = 0.3 and 0.8 given by the Bethe 
ansatz method, respectively. Here we note that the vertical line in the figure indicates the 
eigenvalue E^ divided by U"^ + 4. The eigenstates in a twofold permanent degeneracy in 
the energy spectral flows also permanently degenerate in the spectral flows of I3. Figure 
and Table ini tell us that the eigenvalues E^ of the eigenstates indexed by {— |,— 1,|,|;0;} 
and {— ^,^;;0} never have degeneracies in < n < 1. Hence they belong to the different 
eigenspaces of I^, and the energy level crossing at u = 0.5 is "accounted for" by the third 
dynamical symmetry. Thus all the common eigenspaces of {H, 12,13} with the set of quantum 
numbers {P = 0, S = Sz = r] = r]z = 1} are one-dimensional. 

We give some comments on the results. The origin of the twofold permanent degeneracies 
in Figure \^ has not been clarified until we obtain the correspondence between the energy 
spectral flows and the Bethe states. Indeed we have found that they are due to the reflection 
symmetry. At first sight, there is no crossing at the same value of u in the spectral flows in 
Figure |1J 15] and El But we must not immediately conclude that the energy level crossing in 
Figure m is due to the dynamical symmetries since the spectral flows of I2 and I3 also include 
permanent degeneracies. Actually the eigenstates that have an energy level crossing in Figure 
permanently degenerate in the spectral flows of I2 in Figure [5J Thus the correspondence 
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Figure 4: Spectral flows E from direct diagonalization for L = 6, A?^ = 4, M = 1, P = and 
5 = 77 = 1. 
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Table 3: E obtained by the Bethe ansatz method. 
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Figure 5: Spectral flows E2 from direct diagonalization for L = 6, A?^ = 4, M = 1, P = and 
5 = 77 = 1. 
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Table 4: E2 obtained by the Bethe ansatz method. 
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Figure 6: Spectral flows from direct diagonalization for L = 6, A?^ = 4, M = 1, P = and 
5 = 77 = 1. 
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Table 5: £^3 obtained by the Bethe ansatz method. 
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between the spectral flows and the Bethe states is crucial in the discussion above. 



4 Concluding remarks 

We have studied degeneracies in the energy spectrum of the one-dimensional Hubbard model 
with one down-spin on benzene. The energy spectral flows have been obtained through the two 
approaches: direct diagonalization of the Hamiltonian matrix and the Bethe ansatz method. 
As is noted in IS], the former approach does not always determine the energy spectral flows 
since, if two of eigenvalues approach in numerical data, there are alternative spectral flows 
that may be drawn from the data. To employ the latter, we have presented a procedure 
for giving numerical solutions of the Lieb-Wu equations for UL > 8. Combining the two 
approaches, we have shown that some energy spectral flows have crossings which can not be 
understood by the [/-independent symmetries. We have investigated the spectral flows of the 
higher conserved operators I2 and in the same way and have found that the degenerate 
eigenstates are classified by the dynamical symmetries. 

As we have indicated, the transfer matrix approach systematically provides higher con- 
served operators. Although it is not so easy to see their explicit form, their eigenvalues can be 
written as functions of the solutions for the Lieb-Wu equations (|2.5() by using the eigenvalues 
of the transfer matrix [271 112j . First five of them are given by 

N 

En = y^enjkj), 

i=l 

ei{k) = —2 cos k — -U, 

e2{k) = -2sin(2A;) - 2^7 sin A:, 

3 1 

esik) = 2cos(3A;) -h2cos/!: + C/(3cos(2A;) - l) + -t/^ cos k U , 

4 16 

64 (A;) = 2sin(4A;) + -sin{2k) + uUsm{3k) - - sink) + 2U^sin{2k), 
3 V 3 

65 (A;) = —2 sin(5A;) sin(3A;) sin(3A:) — u(^5 cos k -\ — 

- cos(3A:) - ^ cos A.) - U' cos(2fc) + 1) + If/^ cos - ^U'. 

Our numerical solutions for the Lieb-Wu equations (|2.5() immediately gives their values for 
LU > 8. We have also verified that these eigenvalues never degenerate simultaneously. 

Our studies on the spectral flows of conserved operators support several assumptions that 
have been believed for the Hubbard model and other Bethe-ansatz solvable models. One is 
the one-to-one correspondence between solutions of the Lieb-Wu equations and linearly inde- 
pendent eigenstates. We have found that all the common eigenspaces of conserved operators 
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{/„} are one-dimensional in several subspaces characterized by the sets of quantum numbers 
{N, M, P, S, rj}, which shows linear independence of the eigenvectors therein. Another is the 
algebraic independence of the conserved operators {In}- As is easily found, two commutative 
matrices which have the same type of degeneracy are not algebraically independent since one 
of the two matrices can be expressed by a linear combination of powers of another. In our 
situations, the energy spectral flows have crossings at several values of U which do not pro- 
duce degeneracies for the other conserved operators /„. This fact gives a necessary condition 
for the algebraic independence of H and /„. 

We also remark that our discussion on the Lieb-Wu equations in Section 12.31 give an 
evidence for the validity of the counting of states in Takahashi's string hypothesis. It is 
interesting to generalize our procedure to that with two or more down-spins, which enables 
one to analyse the energy level crossings in Heilmann and Lieb's situations jHj. 
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A Numerical solutions of Lieb— Wu equation 

Following the procedure demonstrated in 12.31 we give the numerical solutions {fcj; A} of the 
Lieb-Wu equations (H^J) for L = 6, iV = 4, M = 1, P = and S" = ?? = 1. Here, for the 
solutions indexed by {/*^*^; J; J'} = {—/*•*''; — J; — J'}, that is, {kf, A} = {—ki; —A}, the Lieb- 
Wu equations ()2.5|) decouple as follows; since "^f^iki = 0(mod27r) and A = 0, the real 
wavenumber solutions can be given by {ki, k2,27r — k2,2TT — ki; X = 0} satisfying 

sin = cot (3/ci) , (i = 1, 2), 

and the fc-A-string solutions are {ki, A:2, tt — i^, vr -|- i^; A = 0} satisfying 

U , . , , U sinh(6C) 

sin /cj = — cot (3/cj I, (i = l,2), sinhf = -, — r . 

4 y 'J' ' ^' ^4 cosh(6C) - 1 

We observe below that the solutions indexed by { — |, — |, |, |; 1; } and {— ^, |, |; ; 1} share 
the same A. The similar situation is seen in other pairs of solutions. Indeed they are in 
the relation of complementary solutions by Woynarovich that is, they give 8 (distinct) 
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solutions of the 8th order algebraic equation 

x^{x^ - 2i(A + 1C//4) - 1) - {x^ - 2i(A - 1C//4) - 1) = 0, 

through X = e^^ . 
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1;} 


4.078300853771271 
0.7z717D7Uoo7ooill 


6.124165064590864 
1 .Oo67z7yo7Dzl4z7 


1 . 131437ooo77oooz 


l 2' 


115. 

2' 2' 2' 


-1;} 


4.646457319558159 

U. ioyUzUz4zooo ( ZZO 


5.556008598803976 
z . zU4oo44o34Uoo i 


— 1 . i3i437ooo776ooz 


1 ^ 
l 2' 


113. 
2' 2' 2' 


0;} 


5.031918491061250 
U.obll4zlDZzL)olJoUo 


5.922043144971536 

i.ZOlzDDololloOOY 


1 n-16 
lU 


{/(*) 


J; J'} 




ki 

c 


k2 


A 


I 2' 


i;;o} 




3.938730108910538 
3.141592653589794 


2.344455198269048 
0.6481740825183825 


-1 X 10-15 


l 2' 


§;;-!} 




3.349805335791242 
4.032164140427823 


1.152236997712285 
0.9264247935487532 


-1.131437583776552 


t 2' 






5.130948309467302 
2.251021166751764 


2.933379971388344 
0.9264247935487532 


1.131437583776552 


I 2' 


i;;o} 




5.03191849106125 
3.141592653589793 


1.251266816118337 
0.6481740825183825 


10-16 


{-§> 






4.833550063047541 
3.763829876294069 


0.2051607987234935 
0.7358326337961154 


-0.7390790775734128 


{-i 


§;;!} 




6.078024508456093 
2.519355430885517 


1.449635244132045 
0.7358326337961154 


0.7390790775734128 








5.922043144971536 
3.141592653589793 


0.3611421622080508 
0.6481740825183825 


10-16 



24 



• n = 0.8 



{/«;J;J'} 




h 


k2 












ki 


A 


{- 


5 3 3 5. 

2' 2' 2' 2' 


0;} 


3.709764868570788 
1.489396519852559 


4.793788787327027 
2.573420438608797 


10-16 


{- 


5 3 15. 
2' 2' 2' 2' 


1;} 


3.950200892165781 
0.7483113471533719 


5.005557244521997 
2.862301130518022 


3.876446955733767 


{- 


5 13 5. 
2' 2' 2' 2' 


-1;} 


3.420884176661565 
1.277628062657589 


5.534873960026214 
2.332984415013805 


-3.876446955733767 


{- 


5 115. 

2' 2' 2' 2' 


0;} 


3.709764868570789 
0.4849278293118341 


5.798257477867752 
2.573420438608798 


10-16 


{- 


5 113. 
2' 2' 2' 2' 


1;} 


3.96270066619563 
0.7650646846850695 


6.03997184302193 
1.798633420456542 


4.232098651896293 


{- 


3 115. 

2' 2' 2' 2' 


-1;} 


4.484551886723044 
0.2432134641576561 


5.518120622494517 
2.320484640983956 


-4.232098651896293 


{- 


3 113. 

2' 2' 2' 2' 


0;} 


4.793788787327027 
0.4849278293118341 


5.798257477867752 
1.489396519852559 


10-16 


{/«;J;J'} 




ki 


k2 










c 




A 


{- 


5 5..Q1 

2 ' 2 ' ' J 




3.70976486857079 
3.141592653589792 


2.573420438608798 
2.094719300551342 


5 X 10-^5 


{- 


5 3.._i-. 

2 ' 2 ' ' J 




3.404556288703187 
3.947809577880995 


1.266195169893996 
2.454678139195706 


-4.232098651896293 


{- 


3 5..;l| 

2 ' 2 ' ' J 




5.016990137285591 
2.335375729298592 


2.8786290184764 
2.454678139195706 


4.232098651896293 


{- 


3 3..QI 

2 ' 2 ' ' J 




4.793788787327028 
3.141592653589793 


1.48939651985256 
2.094719300551342 


2 X 10-^^ 


{- 


3 i.._n 

2 ' 2 ' ' J 




4.503384094618568 

3.903260585058044 


0.2564653496245172 

2.410923664087426 


-3.876446955733767 


{- 


1 3..;l| 

2 ' 2 ' ' J 




6.026719957555069 
2.379924722121542 


1.779801212561019 
2.410923664087426 


3.876446955733767 


{- 


1 i.-ol 

2 ' 2 ' ' J 




5.798257477867752 
3.141592653589793 


0.484927829311834 
2.094719300551341 


10-16 
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